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ABSTRACT 

Two habitable planetary states are proposed: an aqua planet like the Earth 
and a land planet that has a small amount of water. Land planets keep liquid 
water under larger solar radiation compared to aqua planets. Water loss may 
change an aqua planet into a land planet, and the planet can remain habitable 
for a longer time than had it stayed an aqua planet. We calculate planetary 
evolution with hydrogen escape for different initial water inventories and different 
distances from the central star. We hnd that there are two conditions to evolve 
an aqua planet into a land planet: the critical amount of water on the surface 
(Mini) consistent with a planet being a land planet, and the critical amount 
of water vapor in the atmosphere (M^v) that dehnes the onset of the runaway 
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greenhouse state. We find that Earth-size aqua planets with initial oceans < 10% 
of the Earth’s can evolve into land planets if = 3 m in precipitable water 
and Mjni = 5 % of the Earth’s ocean mass. Such planets can keep liquid water 
on their surface for another 2 Gyrs. The initial amount of water and M^v are 
shown to be important dividing parameters of the planetary evolution path. Our 
results indicate that massive hydrogen escape could give a fresh start as another 
kind of habitable planet rather than the end of its habitability. 

Subject headings: planets and satellites: terrestrial planets - planets and satel¬ 
lites: atmospheres - planets and satellites: oceans 


1. Introduction 


The advancement of observation technologies has facilitated an increase in the number 
of extrasolar planets detected. Some of these planets are considered to be terrestrial or rock y 
planets based on the relationship between their mass and radius fe.g.. iRatalha et ahll2013l) . 
Although it is not clear how much water exists on these planets, the marked variations in 
their bulk densities imply that considerable variation exists in the amount of water on these 
extrasolar t errestrial planets. T he habitability of extrasolar planets has been actively dis¬ 


cussed fe.g.. lBorucki et al.ll2013l) . In particular, the stability of liquid water on the planetary 


surface has been extensively examined as liquid water is considered to be essential for the 
origin and evolution of life. 

It is typical to assume that habitable planets, like the Earth, have liquid water oceans. 
Discussions of planetary habitability and the duration of planetary habitability (known as 
continuous habitability) have focused on the stability and evolution of such oceans. If the 
insolation that a planet covered with an ocean receives from a central star is not too high, 
then the upper atmosphere of such a planet would have a very low mixing ratio of water vapor 
due to the presence of a cold trap. Under these conditions hydrogen escape is rather slow, 
and partly dependent on the presence of other gases containing hydrogen (e.g. methane). 
As the luminosity of the central star increases with time, the planet gets warmer and the 
mixing ratio of water vapor in the upper atmosphere increases. When the central star gets 
bright enough, the cold trap disappears and the mixing ratio of water vapor in the upper 
atmosphere becomes high e nough to cause rapi d water loss. This planetary state is called 
the moist greenhouse state (iKasting et al.l 1199311 . Somewhat later, the star becomes bright 
enough to trigger the runaway greenhouse effect: because water vapor is a strong greenhouse 
gas, a posi tive feedback works until all liquid water evaporates in the runaway greenh ouse 


state fe.g.. lAbe fc Matsuilll988l: iKastindllQSSi: IZahnle et al.lll988l: iNakaiima et al.lll992[l . 
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The habitable zone (HZ) has been dehned as th e region aronnd the central star where 


liqnid water is stable on the planetary snrface (e.g., iKasting et alJll993l: iKopparapu et al. 
20131) . The HZs for various pl anetary masses around central stars with various spectral 


types are also investigated (e.g.. IKopparapu et al.ll2014j) . Two concepts have been proposed 


to describe the inner edge of the habitable zone. One is to dehne it as the inner edge 
of the instantaneous habitable zone; that it is, it does not matter how long the planet 
remains habitable. The other dehnes it as the inner edge of the continuously habitable zone 
(GHZ), by which it is meant that a given planet has been in the habitable zone from the 
time it was formed until today. The onset of the runaway greenhouse effect determines the 
instantaneous habitable zone. When considering the long-term stability of liquid water on 
the planetary surface, the water loss from the planet is an important consideration. Under 
moist greenhouse state, hydrogen escapes rapidly. If the escape of hydrogen occurs rapid 
enough, then the inne r edge of the GHZ may be characterized b y the onset of the moist 
greenhouse state fe.g.. iKasting et al.l 1 19931: IKopparapu et al.ll2013l) . 


It has recently been shown that the climates of planets wit h small or large amounts of 
water will differ markedly f rom e ach other ( Abe et ^boOS . 2011 ). Using a general circulation 
model (GGM), lAbe et al.l (120111) examined the habitable zones of land planets. They found 
that the local balance between the precipitation and evaporation of water is controlled by the 
distribution of liquid water on the planetary surface. Liquid water typically accumulates in 
the cooler regions of a planet, that is, high-latitude regions and low-latitude regions typically 
become dry. Given t hat there is no upper limit on planetary radiation in dry areas (e.g., 
Nakaiima et al.lll992[) . a land planet can maintain liquid water on its surface at much larger 
incident of stellar radiation than is possible on aqua planets. In addition, the mixing ratio 
of water vapor in the upper atmosphere is very low on land planets, which means that the 
hydrogen escape flux should also be very low. However, the liquid water on a land planet 
will evaporate if the incident stellar radiation is high enough. This radiation is calle d the 
threshold of the runaway greenhouse for a land planet. According to lAbe et al.l (l201ll) . this 
is 415 W/m^. The inner edge of the habitable zone for a land planet is determined by this 
threshold, which corresponds to 0.77 AU in our solar system. For a.qua p lanets, the onset 
of the runaway greenhouse effect is about 375 W/m^ ( Leconte et ahl 12013 ). Therefore, the 
inner edge of the habitable zone for a land planet is much closer to the central star than that 
for an aqua planet. 


The only difference between a land planet and an aqua planet is the amount of water 
that the planets have. If an aqua planet evolves to a land planet by rapid water loss before 
all of the water on its surface evaporates, then the lifetime of liquid water on its surface can 
be extended because land planets are highly resistant both to water loss and to complete 
evaporation. In this way, rapid water loss might lead to maintaining, rather than ending. 

































4 


the habitable world. 


When lAbe et al.l (1201 ll) discussed the evolution from an aqua planet to a land planet, 
they assumed that the aqua planet can evolve to the land planet when the planetary surface 
becomes almost completely dry before the onset of the runaway greenhouse made the planet 
uninhabitable. They estimated how much water could escape from a planet in Earth’s or¬ 
bit using models that assumed both constant (0.3) and variable albedos. For the constant 
albedo scenario, an aqua planet with Earth’s ocean mass did not lose all of its water before 
the onset of the runaway greenhouse, because the duration of the moist greenhouse state 
was too short. Such planets lapsed into the runaway greenhouse state and became unin¬ 
habitable. Conversely, in the variable albedo scenario, the increase of the albedo due to 
Rayleigh scattering prolonged the duration of the moist greenhouse state. In this case, the 
aqua planet with Earth’s ocean mass lost almost all of its surface water before the onset 
of the runaway greenhouse state, and then evolved into a land planet. Here, we investigate 
the evolution from the aqua planet mode to the land planet mode invoking more realistic 
transition conditions for the amount of water vapor in the atmosphere and the amount of 
water on the planetary surface. 

In the evolution from the aqua planet to the land planet, there is a race between rapid 
water loss and the increase in the luminosity from the central star that triggers the runaway 
greenhouse effect. Whether or not an aqua planet safely negotiates the evolutionary path 
from the aqua p l anet mode to the land planet mode is determined by the outcome of this 


race. 


Abe et al.l fl201ll) only considered the evolution from aqua planet mode to land planet 


mode for planets that have one Earth ocean mass in either the Earth’s orbit or Venus’s 
orbit. However, discoveries of many extrasolar planets have revealed the existence of a great 
diversity of planet types, and we consider that the amount of water that planets have should 
reflect this diversity. 


We therefore systematically investigated evolutionary scenarios in which aqua planets 
can evolve into land planets under conditions with different initial amounts of water and 
distances from the central star. We discuss the evolution from aqua planet mode to land 
planet mode and re-evaluate the inner edge of the HZ. In addition, we investigate the effect 
of varying the transition conditions that determine whether an aqua planet can evolve into 
a land planet. 
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2. Evolution model 


In this paper, we consider three planetary states (a water planet state, a steam planet 
state, and a dry planet state), and two modes for a water planet state (an aqua planet mode 
and a land planet mode). Here, we describe the main characteristics of these states and 
modes in advance. 


In a water planet state, liquid water is stably present on the planetary surface. In most 
previ ous studies, planets i n the aqua planet mode are considered to be potentially habitable 
(e.g., iKasting et al.lll993l: iKopparapu et al.l 120131) . In these studies, it has been implicitly 
assumed that planets have a large amount of water (e.g., the present Earth’s oceans). On 
the other hand, if planets have a s mall amount of water, then they are considered to be in 
land planet mode flAbe et al.ll201l[) . 


Planets in a steam planet state have water vapor in their atmosphere but no liquid water 
on their surface. If planets with liquid water receive the radiation exceeding the threshold 
for the runaway greenhouse effect, then all of the liquid water on their surface will evaporate 
and these planets will evolve into a steam planet state. 


In a dry planet state, planets have no liquid water on their surface and no water vapor 
in their atmosphere. Dry planets are not the same as land planets, as the latter have a 
small amount of liquid water on their surface. If water vapor escapes from a steam planet, 
then it becomes a dry planet and both steam planets and dry planets are considered to be 
uninhabitable. 


Here we simulate the changes of the amount of water on the planetary surface and in 
the atmosphere using a numerical model. In section 2.1, we describe the basic processes that 
cause changes in the amount of atmospheric and surface water. The transition conditions 
from aqua planet mode to land planet mode are described in Section 2.2. In Section 2.3, 
we describe the hydrodynamic escape of hydrogen as a rapid water loss. The escape flux 
of hydrogen molecules depends on the mixing ratio of hydrogen molecules in the upper 
atmosphere and the evolution of the luminosity and EUV flux of the central star, which are 
described in sections 2.4 and 2.5, respectively. 


2.1. The processes that change the amount of water 

The initial amount of water on terrestrial planets just after their formation is not known. 
We treat the initial amount of water as a parameter, i.e., we consider planets with various 
initial amounts of water. The amount of water on the surface of a planet would change 
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through the interaction with the planetary interior through processes such as subduction 
of hydrous minerals, outgassing of water at mid-ocean ridges and arcs, and so on. For the 
present Earth, the ingass i ng flu x of water is estimated to be 5 x 10^^ — 1.1 x lO^^mol/Gyr 
( Rea fc Rufi 19961: Javov 1998h. and the degassing flux of wate r is estimated to be 2.2 x 
1022_i 1 X10^^ mol/Gyr flBounama et al.ll200ll:lHilton et al.ll2002l) . Ingassing and outgassing 
fluxes in water appear to be in equilibrium on the present Earth. The difference between the 
ingassing and outgassing water fluxes is the order of 10^^ mol/Gyr. By comparison, once the 
rapid escape of water into space starts as the central star evolves, the escape flux of water 
vapor exceeds 10^^mol/Gyr(see Sections 3.1 and 3.2). For simplicity, we assume that the 
change in the total amount of water on a planet is attributed solely to the loss of water into 
the space. 

We assume that the amount of water decreases via the hydrodynamic escape of hydrogen 
molecules. The change in the amount of water is given by 


dt 


( 1 ) 


where MhjO is the total amount of water and 0esc is the escape flux of molecular hydrogen. 
The escape process of hydrogen is described in Section 2.3. 

The amount of water vapor in the atmosphere depends on the atmosphere model used. 
For the aqua planet and land planet modes, the amount of water vapor is esti mated by using 
the 1- D, cloud-free, radiative-convective atmospheri c mod el developed by I Abe fc Matsui 
(1l988l) and using the results obtained by I Abe et al.l (120111) . respectively. The amount of 
water on the surface can then be estimated by subtracting the amount of water in the 
atmosphere from the total amount of water on the planet. An aqua planet is considered to 
evolve into a land planet if the amounts of water in both reservoirs meet specihc transition 
conditions, described in detail in Section 2.2. 


Planetary albedo, which plays an important role in climate, depends on the properties 
of the planetary surface and atmosphere. In our numerical simulation, we hxed the albedo 
of the planetary surface at 0.3, which is typical of deserts. When the amount of water 
vapor increases due to warming caused by an increase in the luminosity of the central star, 
the atmosphere gets denser and light from the central star is scattered more efficiently. 
Therefore, we also consider the scattering albedo for an aqua planet. Since the scattering 
albedo is negligible for a land planet (if the effect of dust is not considered), we consider that 
the planetary albedo for a land planet to be 0.3. 

Here, we assume that the planets have zero obliquity and eccentricity, with a mass 
and a radius equivalent to those of the Earth. In the planetary atmosphere, 1 bar of the 
present Earth’s atmosphere is considered as the background atmosphere. Except for water 
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vapor, no changes in atmospheric components were considered in this stndy. It is possible 
that the carbon cycle exists on Earth-like planets. The carbon cycle in the presence of 
liquid w ater provides a neg ative feedback on surface temperature that maintains a stable 
climate flWalker et al.lIlQSlIl . However, in the case of a central star with high luminosity, the 
feedback breaks down because there is too little CO 2 in the atmosphere to affect changes in 
surface temperatures owing to efficient chemical erosion (rapid hxing of CO 2 as carbonate 
rock). Therefore, the evolution of the atmosphere is driven almost entirely by changes in the 
amount of water in the planetary atmosphere. 


2.2. Transition conditions 


Abe et al.l (120111) assumed that an aqua planet with the present Earth’s ocean mass could 
evolve to land planet mode just before all of its water is lost, and before there was a rapid 
increase in surface temperatures due to the runaway greenhouse effect. Under such transition 
conditions, it is possible that a planet with a large amount of water in its atmosphere can 
evolve from an aqua planet to a land planet. However, if the planet has a large amount 
of water vapor in its atmosphere before the loss of water, then it may maintain a heated 
state during atmospheric escape. Under these conditions, such a planet may not evolve into 
a land planet, even if it loses most of the water on the surface. It is therefore necessary 
to consider the amount of water vapor in the planetary atmosphere. We therefore imposed 
several conditions on amounts of water on the planetary surface and amounts of water vapor 
in its atmosphere for the transition from an aqua planet to land planet mode. 

However, one problem is that we do not know how much water vapor there is in an aqua 
planet’s atmosphere just before the threshold for the runaway greenhouse effect is reached. 
In order to evolve to land planet mode, there should be an upper limit to the amount of water 
in the atmosphere; we treat this upper limit to the amount of water as a parameter, which 
we denote as Mcv and refer to a s ’’the critical amount of vapor for the runaway greenhouse 
e ffect”. From previou s studies flAbe fc MatsuilIl988l: iKastindIl988l: iKopparapu et ahl 12013 


Goldblatt et al.ll2013l) . it appears that when Mc^, is expressed as a column it is equivalent to 
1 m to 10 m of precipitable water. In this paper, to illustrate typical results, an Mcv of 3 m 
was used for our standard model. 

Another important parameter is the maximum amount of water that a land planet can 
have and still behave as a l and planet. We consider the distribution of surface water from the 
results of I Abe et al.l (120111) . Surface water on a land planet is typically located at latitudes 
above 60 degrees, which corresponds to 10 % of the entire planetary surface. According to 
percolation theory, for any given area, only half of that area is required to form a connected 
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area. The maximum unconnected oceanic area is thus half of the area above 60 degrees 
latitude, which corresponds to 5% of the entire planetary surface. If planetary topography 
is ignored, then the typical amount of water on a land planet is 5 % of the present Earth’s 
ocean mass. If a planet has more water, that water will extend equatorward of 60 degrees 
latitude. Therefore, we assume that 5 % of the present Earth’s ocean mass on the planetary 
surface approximates the transition conditions that are required for the maximum amount 
of water required for an aqua planet to evolve into a land planet. We refer to this condition 
“the maximum liquid water mass for a land planet mode”, which we denote as Mmi. If the 
amounts of water on the planetary surface and in the atmosphere meet these conditions, then 
the aqua planet can evolve into a land planet. We treat Mmi = 0.05 Mqc as our standard 
case, but we also treat Mmi as a parameter that varies from 0.01 to 0.1 Mqc, where Mqc is 
the amount of the present Earth’s ocean (8.4 x 10^^ moles). In Section 3.5, we discuss the 
evolution of an aqua planet under several transition conditions. 


2.3. Water loss 


Escape is an important process in atmospheric evolution. In this study, we consider the 
hydrodynamic escape of hydrogen molecules because it is the fastest process of hydrogen 
escape and has the most potential to decrease the water reservoir of a planet. We will 
address two major bottlenecks to hydrogen molecular escape from the planetary atmosphere. 
The hrst is the energy source required to drive the escape of hydrogen molecules. The 
main source of energy for escape is extreme ultraviolet radiation (EUV), usually dehned 
by A < 100 nm from the central star, which is directly absorbed by hydrogen molecules 
and atoms. When the escape flux is controlled by the EUV flux, the escape is said to 
be “energy-limited”. The other major bottleneck is the diffusion of hydrogen-containing 
molecules through the background atmosphere to reach the altitudes where escape occurs. 
When escape is controlled by the diffusion of molecules through the upper atmosphere, then 
the escape is referred to as “diffusion-limited”. We calculate escape fluxes in both modes 
and adopted the smaller one as the actual escape flux. Other possible bottlenecks, such 
as the photon-limited photochemistry that converts H 2 O to H 2 , magnetospheric drag, and 
radiative cooling by ions (e.g., Hg ) and molecules embedded in the ontflowing wind, are less 
well established and will be neglected here. 


The escape flux of the diffusion-limited escape mod e depends on the total mixing ratio of 
hydrogen-containing species in the upper atmosphere fjWalkerl 119771) . The diffusion-limited 
escape flux ) of hydrogen molecules is given by 

b{m^ - m,)g 


— /t(H2 


kT. 


str 


( 2 ) 
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where /t(H 2 ) is the total mixing ratio of hydrogen molecules in all forms at homopause 
(/t(H 2 ) = /(H 2 O) + /(H 2 ) + 2 /(CH 4 ) + ■•■), Tsti is the temperature in the stratosphere 
(^str = 200 K), b is the binary diffusion coefficient between the background atmosphere and 
H 2 ( fo = 1.9 X 10^®(Tstr/300 cm“^s“^ for H 2 in air), g is the acceleration due to gravity, 

k is Boltzmanns constant, rua and rrii are the average molecular masses of the background 
atmosphere and the molecular mass of hydrogen, respectively. We assumed that the major 
carrier of hydrogen is water vapor. Therefore, we adopt /t(H 2 ) = /(H 2 O), where /(H 2 O) is 
the mixing ratio of water vapor in the upper atmosphere. 


An escape flux is limited by EUV flux if the planetary upp er atmosphere is hyd rogen-rich 
atmosphere. The energy-limited escape flux {(j)e ) is given by IWatson et ahl fjl98ll) as 


e*S'Euv’"p 

rriiGMp 


( 3 ) 


where e is the escape efficiency, S'euv is the EUV radiation flux, G is the gra vitational con 


stant, Tp is the radius of the planet, and Mp is the planetary mass. According to lWatson et ah 


fjl98ll ). the escape efficiency ranges from 0.15 to 0.3; we adopted 0.15 as the escape efficiency. 
Since we used Earth-sized planets in this study, we used the Earth’s radius and mass for rp 
and Mp, respectively. 

Atmospheric H 2 O molec ules are dissociat ed by UV radiation from the central star (via 
H 2 O -|- hi/ ^ H -|- OH) (e.g., iBriukmaunIIl969h . Hydrogen molecules are transported to the 
upper atmosphere, where EUV radiation reaches, and are heated by su ch radiation. The 
resultant hydrogen molecules escape to space(e.g., iKasting et al.Nl984i) . When hydrogen 
escapes, the oxygen left behind can build up in the planetary atmosphere. Because we 
assume diffusion-limited escape, oxygen does not escape. If an amount of w ater equivalen t 
to the Earth’s ocean mass escapes, approximately 240 bars of O 2 will remain (jKastindll997l) . 
In this study, we assume that all of this oxygen is absorbed by the crust of the planet. 


2.4. The mixing ratio of water vapor in the upper atmosphere 


In order to estimate the escape flux of hydrogen molecules in the diffusion-limited escape 
mode, we need the mixing ratio of water vapor in the upper atmosphere. For the aqua planet 
mo de, we estimat e the m ixing ratio of water vapor in the upper atmosphere as described 
by lAbe fc Matsuil fjl988l) . Figure [T] shows the mixing ratio of water vapor in the upper 
atmosphere as a function of the incident stellar radiation. It is very low until the cold trap 
region disappears. The escape flux of hydrogen molecules in the diffusion-limited escape 
mode is thus expected to be very low due to the very low mixing ratio of water vapor under 
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these conditions. If an aqua planet receives incident solar flux from the central star that is 
sufficiently strong to eliminate the cold trap, the upper atmosphere will become humid and 
the escape of hydrogen molecules will be rapid. 


For the land planet mode, we estimate the mixing ratio of water vapor in the upper 
atmosphere using the results of lAbe et ahl (120111) . The mixing ratio for a land planet also 
increases with increasing luminosity of the central star. In a case where the insolation from 
the central star is above the threshold for the runaway greenhouse effect for a land planet, 
the upper atmosphere becomes humid and the escape of hydrogen is rapid (FigJ2]). 


2.5. Stellar evolution 


We consider the evolution of total luminosity and the EUV flux of a star because the 
escape flux of hydrogen molecules depends strongly on both. The total luminosity of a star 
increases over time. We treated the central star as a so lar type star. For the luminosity for 
a G-type star, we use the following expression given bv iGoughI fjl98ll) : 

-1 


Sit) = 


2 

'+5 


t 

1 - 

f© 




0 


(4) 


where S{t) is the incident stellar flux for a G-type star at lAU as a function of time, t is 
the time in Gyr, t© is the age of the Sun ( t© = 4.5 Gyr ), and S'© is the solar constant ( 
Sq = 1366 W/m^ ). 

Since the EUV flu x for young sta r s is st ronger than that for old stars, we employed the 
following expression of iLammer et ahl (1200911 for a G-type star: 


-^EUV — 


(5) 


0.375Lof°-^2^ (t < 0.6Gyr) 

0.19Lor^®9 {t > 0.6Gyr) 

where Teuv is the luminosity of EUV in Watt s“^, and Lq = 10^^'^^W. The relation between 
the luminosity and the flux is given by S'euv = -hEuv/4vra^, where a is a distance from the 
central star to the planet in meters. 


3. Results 

3.1. Effect of the initial amount of water 

In this section, we examine typical results obtained by our numerical simulations de¬ 
scribing the planetary evolution. Figure |3] shows the evolution of a planet with the present 
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Earth’s ocean at the beginning (Mini = 1-0 Mqc) at a = 0.80 AU, where M^i is the initial 
amount of water. 

As long as the luminosity of the central star is low (< 3Gyr), the mixing ratio of 
water vapor in the upper atmosphere remains low (Figure [Sti). Consequently, the escape 
flux of water vapor in the diffusion-limited escape mode also remains extremely low (Figure 
[3b) and the amount of water on the surface hardly changes during this period (Figure [3b). 
The atmosphere gets wetter as the luminosity of the central star increases over time, and 
after 3 Gyr, the mixing ratio of water vapor in the upper atmosphere increases rapidly. The 
increase in the mixing ratio in the upper atmosphere enhances the escape flux of water 
vapor in the diffusion-limited escape mode and decreases the amount of water on the planet 
(Figure [3b)- In this calculation, this planet does not satisfy the transition conditions from 
aqua planet mode to land planet mode because the amount of surface water exceeds the 
maximum water mass for land planet mode (0.05 Mqc) when the average column mass of 
atmospheric water vapor reaches Mcv (i.e. a depth of 3 m), which is the critical amount for 
the runaway greenhouse effect. Therefore, the planet remains in aqua planet mode until the 
onset of the runaway greenhouse effect, eventually evolving into the uninhabitable steam 
planet state. Over a period of approximately 0.1 Gyr, all of the water escapes into space and 
the planet finally evolves into a dry planet state. 

Figure [H shows the evolution of a planet with Mini = 0.1 Mqc and a = 0.8 AU. The 
difference between the cases shown in Figure [3] and Figure [H is the initial amount of water. 
The overall evolution of water is very similar to that shown in Figure [3] before t = 3.19 Gyr. 
However, when the amount of surface water decreases to Mmi = 0.05 Mqc (i.e. the maximum 
water mass for the land planet mode), the average column mass of the atmospheric water 
vapor is less than the 3 m threshold for the runaway greenhouse effect in hgure [Hi. This 
planet therefore satishes the requirements for a transition from aqua planet mode to land 
planet mode and can evolve into land planet mode. Once in land planet mode, the planet’s 
atmosphere becomes dry (Figure [3b) and escape shuts off (Figure [3b)- This planet can 
keep liquid water on its surface until the central star has become sufficiently bright enough 
to exceed the runaway threshold of a land planet (Figure [3b)- When the planet receives 
insolation above this limit, the planet evolves into the steam planet state and then the dry 
planet state once it has lost its remaining water into space over about 0.17 Gyr. 

Whether or not an aqua planet can evolve to a land planet thus depends on the initial 
amount of water. Qualitatively, planets with small initial amounts of water can more easily 
evolve from the aqua planet mode to the land planet mode. 
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3.2. Effect of the distance from the central star 

Here we examine the effect of the distance from the central star on planetary evolution. 
Figure [5] shows the evolution of a planet with Mini = 0.1 Mqc and a = 1.0 AU. In Figure [S]!, 
the amount of surface water remaining on the planet exceeds M^i and the planet experiences 
the runaway greenhouse effect. This planet evolves into the steam planet state without 
passing through a land planet state. On the other hand, we have seen that a planet with the 
same initial amount of water at 0.8 AU can evolve from aqua planet mode to land planet mode 
(Figure Hj). This evolution occurs because the planet nearer the central star evolves through 
the moist greenhouse stage when the star is younger, and the EUV flux and hydrogen escape 
flux are higher. Any planet near the central star will lose more water during the period of 
rapid escape of water vapor. 

The timing of the onset and duration of rapid escape of water vapor are the two keys 
that determine whether or not an aqua planet will evolve into a land planet. In the case of a 
planet that is located far away from its central star, the onset of rapid escape of water vapor 
occurs later when the star is older and its EUV flux is lower. It therefore takes more time 
for a more distant planet to lose its surface water; phrased another way, the more distant 
planet loses less water between the onset of the moist greenhouse effect and the onset of 
the runaway greenhouse effect for an aqua planet. Moreover, because the rate of stellar 
luminosity increases with time (see (jl])), more distant planets have less time to shed excess 
water before the onset of runaway greenhouse effect. As a result, it is harder for a distant 
planet to evolve into a land planet, and easier to evolve directly into a steam planet. 


3.3. Limit of the initial amount of water to evolve to a land planet 

In the previous section, we showed several typical numerical results obtained for plane¬ 
tary evolution. We found that the initial amount of water (Min,) and the distance from the 
central star (a) are important in determining whether or not an aqua planet can evolve into 
a land planet. 

Figure [H] shows the maximum initial amount of water that is required for the transition 
from the aqua planet mode to the land planet mode as a function of the distance from the 
central star. In this hgure, planets below the curve can evolve from aqua planet mode to 
land planet mode. Aqua planets close to a central star with small initial amounts of water 
can evolve into land planets, but aqua planets that are too far from a central star, or that 
are born with large amounts of water, evolve into steam planets without passing through 
the land planet stage. 
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Planets with less than the maximum liquid water mass for land planets (Mini < Mmi) 
remain in land planet mode from the beginning. In cases of planets with small orbits (a < 
0.72 AU), all of the water evaporates and such planets are in a steam planet state from the 
beginning. 


3.4. Inner edge of the HZ 


Here we discuss the inner edge of the HZ of solar-type stars based on our results. Figure 
[7] shows the evolution of planets with Mini = 1-0 Mqc at various distances from the central 
star under the same transition conditions. For planets with Mini = 1-0 Mqc, rapid escape of 
water occurs, but such planets cannot evolve to the land planet mode because they cannot 
lose enough water - they cannot get from 1.0 Mqc to Mmi - before encountering the runaway 
greenhouse threshold for aqua planets. Therefore, the boundary between the aqua planet 
mode and the steam planet state corresponds to t he inner edge of the classical instantaneous 
habitable zone, proposed in iKasting et al.l (119931) (Figure [7]). 


As in Figure [71 Figure [8] shows the time evolution of planets with Mini = 0.1 Mqc. When 
located between 0.72 and 0.82AU, such planets can evolve from aqua planet mode to land 
planet mode. Planets in the land planet mode remain habitable for an additional 2 Gyrs 
before evolving into the steam planet state and then to the dry planet state due to increased 
luminosity of the central star. 


In both Figure [7] and Figure [HI aqua planets closer than about 0.7 AU from the central 
star are in the steam planet state from the beginning. Such planets lose all of their water 
rapidly (over a few tens of million years) because t he EUV flux at that time is very strong; 
such planets are called Type H terrestrial planets ([Hamano et al.[[2013[) . 


Figure [9] shows the inner edge of the GHZ, which is dehned here as the region where 
a planet has liquid water on its surface for more than 4.6 Gyr. We divided the GHZ into 
two types according to the initial amount of water. One type of GHZ refers to those regions 
on aqua planets that maintain wa ter on their su r face f or more than 4.6 Gyr; this type cor¬ 
responds to the GHZ proposed in [Kasting et al.[ fjl993[) and is referred to here as a type-I 
GHZ. The other type of GHZ refers to regions on land planets that maintain water on their 
surface for a total of more than 4.6 Gyr; this type is referred to as a type-H GHZ. 


In our model, planets with Mini > 0.22 Mqc cannot evolve to land planets: such planets 
have a type-I GHZ. On the other hand, planets with M^i < 0.22 Mqc can evolve into land 
planets. However, planets with 0.19 Mqc < M^i < 0.22 Mqc cannot maintain liquid water on 
their surface for 4.6 Gyrs because the insolation on such planets exceeds the threshold of the 
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runaway greenhouse effect. Planets with Mini < 0.19 Mqc have a type-II CHZ. 

Figure [9] shows that planets that begin with 0.09 Mqc < M^i < 0.19 Mqc have both 
type-I and type-II CHZs, but these two CHZs are not contiguous. However, planets with 
0.05 Mqc < Mini < 0.09 Mqc also have both CHZs and they are contiguous. Planets with an 
initial water less than the maximum water mass for a land planet (Mmi = 0.05 Mqc) on their 
surface are in land planet mode from the beginning. 


4. Discussion 

4.1. Effect on the transition condition 

In previous sections, we showed results when Mcv, the critical water vapor amount for 
the runaway greenhouse, was taken as a global average of 3 m water depth equivalent and 
Mmi, the maximum liquid water mass for a land planet mode, was taken to be 0.05 Mqc. In 
this section, we investigate how the evolution from aqua planet mode to land planet mode 
changes when the transition conditions are altered. 

Figure fTOl shows the maximum M^i required for a planet to evolve into a land planet 
for different values of Mcv. In these calculations, M^\ is fixed at 0.05 Mqc. If Mcv exceeds 
3 m, then the maximum M^i increases because the onset of the runaway greenhouse effect is 
delayed. For planets closer to the central star, the maximum Mini is larger because the escape 
flux of water is larger before the onset of runaway greenhouse effect. However, if planets are 
too close to the central star, the maximum M^i decreases and has a peak (see Figure fTOj) . 
The amount of water that escapes is determined by the duration of the period of a rapid 
escape until the onset of the runaway greenhouse effect and the escape flux of water during 
this period. For planets that are very close to the central star, runaway greenhouse effect 
occurs early and duration of the period of rapid escape is short. Although the escape flux 
in the energy-limited escape mode is large, the rapid escape of water occurs in the diffusion- 
limited escape mode for these planets. Therefore, the total amount of escaped water hits a 
peak, which is shown by the maximum Mini value in Figure [TOl An aqua planet cannot evolve 
into a land planet when M^ is 1 m (global average). In our model, the mixing ratio of water 
vapor in the upper atmosphere reached 10“^ when the column of water in the atmosphere 
is 1.9 m (global average). If M^ is less than 1.9m, then a planet in the aqua planet mode 
cannot lose enough water and it passes directly into the steam planet state. 

Our standard case takes M^i = 0.05 Mqc. This value of M^i is dependent upon the 
planetary topography. Figure fTTl shows the dependence of the maximum Mini on Mmi for 
a planet evolving into a land planet. In this calculation, the value of Mcv is fixed at our 
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standard value of 3 m. A smaller value for Mmi results in smaller maximum Mini, but it is 
possible for aqua planets to evolve into land planets even if M^\ = 0.01 Mqc. However, an 
aqua planet cannot evolve into a land planet when Mcv is less than Im (see Figure [HI]) . 
Therefore, the parameter M^ has more impact on the evolution path to land planets than 
the parameter Mmi. 


4.2. Atmospheric models 


We investigated the evolutionary paths of plane ts using a 1-D, clou d-free, radiative- 
convective atmospheric model developed in the 1980s fjAbe fc Matsuilll988[) . The 1-D atmo- 
spheric models have been improved by using more accurate opacity estim ates for water vapor 
(iKopparapu et al.l 120131: iGoldblatt et al.l 120131: iKopparapu et al.l 120141). The inner edge of 
the GHZ (dehned by moist greenhouse effect) in IKopparapu et al.l f|2013|) is 0.99 AU. In our 
model, the onset of the moist greenhouse effect is 1.39 Sq. Since the onset of the runaway 
greenhouse effect is dependent on the val ue of M^-y, when M , y is se t to 3 m, the onset of the 
runaway greenhouse effect is 1.40 S'©. In IKopparapu et ahl f|2013]) . the onsets of the moist 


and runaway greenhouse effects are 1.015 and 1.06 S'©, respectively. The timing of the onset 
of rapid water loss and the amount of water that escaped should therefore be different be¬ 
tween our atmospheric model and their model. In order to investigate the evolution from an 
aqua planet to a land planet, we need the relationship between the incident solar flux and 
the amount of water vapor in the atmosphere, but it was n ot shown in the pr e vious papers. 
We therefore qualitatively compared our results to those of IKopparapu et al.l f|2013l) . 


We estimated how much water escapes during the period from the onset of the moist 
greenhouse effect to th e onset of the runaway greenhouse effec t using the results from 


Kopparapu et al.l f|2013]) . To estimate the escaped water amount bv IKopparapu et al.l f|2013]) . 
we presume that the mixing ratio of water vapor is 10“^, so that the escaped water amount 
during this period in the diffusion-limited escape mode is approximately 3 x 10^^ moles (ca. 
0.04 Moc). The amount of escaped water under energy-limited escape mode conditions in 
the same period was estimated to be 10^^ moles (ca. 1.19 Mqc). The actual amount of es¬ 
caped water falls between these estimates (i.e., between 3 x 10^^ and 1 x lO^^moles). In our 
model, planets that can evolve into land planets typically lose ab out 0.1 Mnr of wa t er (se e 


Section 13.31). which is co mparable to, or less than, the estimates of IKopparapu et al.l (120131) . 


Kopparapu et al.l (j2013[) indicates that planets at a < 0.84 AU are in the steam planet state. 


Thus, the region where aqua planets can e v olve i nto land planets moves outward (about 0.1 
AU) when the model of IKopparapu et al.l f|2013l) is applied. The type-II GHZ (see Figure 
[9|), in which aqua planets evolve to land planets, also moves outward. 
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4.3. Intensity of EUV flux 


We assessed the effect of EUV ffux on ou r results. The uncertainty associated with th e 
EUV ffux from the central star is large (e.g., IZahnle fc Walken Il982l: iT^ammer et al.ll2009l) . 
Just before an aqua planet evolves into a land planet, water escapes rapidly in the energy- 
limited escape mode (see Figure HJd). If the EUV ffux is larger than that used in our model, 
then the escape ffux increases. However, since the escape ffux is limited by the diffusion 
process, a larger EUV ffux does not markedly affect our results. On the other hand, if the 
EUV ffux is smaller than that used in our model, then the escape ffux decreases. This makes 
it harder for aqua planets to evolve into land planets. 


4.4. Amount of atmospheric CO 2 


An important complication could be the evolution of the amount of CO 2 in atmosphere. 
Since the carbon cycle on a land planet has not been investigated yet, we speculate possible 
effects of carbon cycle on a land planet climate in the following. 

The amount of CO 2 in the atmosphere is governed by the balance between the degassing 
ffux of CO 2 and the removal ffux of CO 2 through chemical weathering. We expect the 
degassing ffux of CO 2 to be smaller on a land planet than on an aqua planet because plate 
tectonics is likely to be less efficient on a water-poor mantle. Chemical weathering may be 
also less efficient on a land planet, because it occurs mainly in wet regions. If the former 
effect exceeds the latter, then the amount of CO 2 in the atmosphere should be smaller on a 
land planet than on an aqua planet. However, any decrease in the degassing ffux would be 
difficult to estimate. 


Therefore, it is not clear whether there would be more or less CO 2 on a land planet 
compared to an aqua planet. If the amount of CO 2 is balanced at a high value, then the 
period of habitability may be shorter than that estimated in our study because the greenhouse 
effect attributed to CO 2 may trigger the runaway greenhouse effect. On the other hand, if 
the CO 2 amount is balanced at a low value, the environment is cooler than our estimate. 
However, such an effect is likely minor, because the insolation is high enough to keep from 
the global freezing without CO 2 greenh ouse effect. It sho uld also be noted that a land planet 
is relatively resistant to global freezing flAbe et al.ll201ll) . Even if ice caps appear, their size 
is regulated by the carbon cycle, because CO 2 removal through chemical weathering would 
be negligible in the ice-covered area. Under such a scenario, the lifetime of a habitable 
environment would likely be similar to that estimated in our study. 
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4.5. The geologic H 2 O cycle 

We have neglected the geologic water cycle by presuming that the fluxes of water into 
and out of the mantle are in rough balance. This rough balance is likely to break down after 
an aqua planet has become a land planet, because the mantle may still be wet and out- 
gassing may still be considerable, but subduction of hydrous minerals is probably restricted 
to locations near the poles. Thus water outgassing would be a potential hazard to the future 
habitability of the land planet over long time scales. But while the aqua planet is losing its 
hydrogen in the moist greenhouse state, it is possible that the imbalance goes in the other 
directions, as the generally warm wet climate of the moist greenhouse would likely promote 
weathering reactions and possibly promote subduction of hydrous minerals. Thus we think 
it probable that the geologic water cycle would aid rather than subvert the transition of a 
suitable planet from aqua to land modes. 


4.6. Other stars 


Stars of different spectral types evolve differently from the Sun. In particular, the 
luminosity evolution of stars of later spectral type (late G dwarfs, K dwarfs, and M dwarfs) 
is slower than in the Sun. These stars are fainter than the Sun and hence the HZ is closer 
to these stars than to the Sun. The ratio of EUV radiation to lumino sity is about the sarn e 
for K stars as for G stars, and the ratio is higher for most M stars flTjammer et al.ll2009l) . 
Hence these stars are at least as effective at driving hydrogen escape as a Sun-like star, and 
so we would expect similar or faster rates of hydrogen escape for an HZ planet around these 
other stars. The slower rate of luminosity evolution means that there is more time available 
for escape to take place before the moist greenhouse becomes a runaway greenhouse, and 
therefore we expect that late type stars are more favorable to a planet making a continuously 
habitable transition from aqua state to land state. We will quantitatively investigate the 
evolution path from the aqua planet mode to the land planet mode for stars of different 
spectral types in the next paper. 


5. Conclusion 

Planets can evolve from aqua planet mode to land planet mode and, in so doing, main¬ 
tain liquid water on their surface for a longer time. Here we assess the various requirements 
of such a transition. We argue that the transition from an aqua planet to a land planet can 
occur if the amount of water on the surface is reduced rapidly by hydrogen escape via the 
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photodissociation of H 2 O vapor, while the amount of water vapor in the atmosphere remains 
sufficiently low to prevent the onset of runaway greenhouse effect; i.e., the planet must lose 
most of its water while in the moist greenhouse state. We considered two factors as being 
central to the evolution of a land planet from an aqua planet: the critical vapor amount for 
the runaway greenhouse, Mcv, and the maximum liquid water mass for a land planet mode, 
Mmi. We assumed that the star in our simulation was similar to the Sun. We examined 
temporal changes in the amount of water on a planet under conditions of different initial 
amounts of water and orbital distance. We illustrated the conditions for the evolution of an 
aqua planet to a land planet and applied this information to re-evaluate the inner edge of 
the HZ. 

The evolutionary paths of aqua planets are mostly determined by the initial amount of 
water. Mini. Planets with large Mini evolve into the runaway greenhouse state before entering 
the land planet mode, while planets with small M^i can evolve into land planets. Planets 
nearer the central star can lose water more quickly because these events take place when 
the star is younger and thus a stronger source of the EUV radiation that drives hydrogen 
escape. These Endings imply that it is generally easier for a land planet to form when it is 
closer to its star than when it is far from its star. 

If Mcv is equivalent to 3 m of precipitable water, then an aqua planet with Mini = 0.1 Mqc 
at 0.8 AU can evolve into a land planet. In addition, such a planet can retain liquid water on 
its surface for an additional 2 Gyr as a land planet. On the other hand, if Mcv is less than 
about 1 m, it is very difficult for any aqua planet to evolve into a land planet. This value 
for Mcv corresponds to the amount of water vapor in the atmosphere that causes the rapid 
water loss of a substantial fraction of the Earth’s ocean mass. By contrast, the dependence of 
the evolutionary path on M^i (the upper bound on the amount of water that defines a land 
planet) is not very large, implying that the value of Mcv is more important as a transition 
condition from an aqua planet to a land planet. 

We describe a new type of continuous habitable zone (GHZ) where aqua planets can 
evolve into land planets and maintain liquid water on their surface for more than 4.6 Gyr(type- 
II GHZ in Figure [H]). Rapid escape of water is not necessarily the end of a planet’s habit¬ 
ability; rather, such water loss can bring about a new and different kind of habitability. 
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Fig. 1.— The mixing ratio of water vapor in the upper atmosphere and scattering albedo 
by water vapor for an aqua planet as a function of the incident solar radiation. The mixing 
ratio of water vapor increases rapidly at around 470 W/m^. In our model, we assume a 
constant albedo for the planetary surface (0.3) and the scattering albedo by water vapor, 
which varies with the amount of water vapor in the atmosphere. 



O 

■D 

0 

n 

CO 

O) 

c 

"l— 

0 
-I—> 
-I—' 

CO 

o 

CO 




22 



350 400 450 500 550 600 650 
Incident solar radiation [ W / ] 


Fig. 2.— The mixing ratio of water vapor in the atmosphere for a land planet as a function 
of the incident solar radiation. If a land planet receives the incident solar radiation exceeding 
the threshold of the runaway greenhouse (585 W/m^) for a land planet, then the mixing ratio 
of water vapor is high. In our model, we consider the constant albedo (0.3) throughout a 
period of time of a land planet because there is little water vapor before the threshold of the 
runaway greenhouse for a land planet. 
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Fig. 3.— Evolution of a planet with Mini = I-OMqc at a = 0.80 AU. (a) The mixing ratio of 
water vapor in the upper atmosphere; (b) escape flux of hydrogen; (c and d) change in the 
amount of water in atmosphere and on the surface are shown as a function of time. Figure 
(d) shows the change in the amount of water around 3.2 Gyr with two transition conditions 
(the critical amount of vapor for the runaway greenhouse effect, Mcv, and the maximum 
liquid water mass for a land planet, Mmi) from an aqua planet to a land planet. Here, M^v 
are set a 3 m in global average for the precipitable water and M^i is set to 0.05 Mqc. In this 
case, the aqua planet evolves directly to a steam planet and then a dry planet; it is never a 
land planet. 
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Fig. 4.— Like Figure El except that the planet has less initial water, Min, = 0.1 Mqc at 
a = 0.80 AU. In this case, the aqua planet can evolve to a land planet and maintains liquid 
water on its surface for an additional 2 Gyrs. 
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Fig. 5.— Like Figure IU but for a planet a little farther from the central star, with Mini = 
0.1 Mac at a = 1.0 AU. Figure (d) shows the change in the amount of water at approximately 
7.7 Gyr. This aqua planet does not pass through a stage as a land planet, but instead evolves 
directly into a steam planet and then a dry planet. 
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Fig. 6.— Initial planetary conditions required for the evolution of a land planet with a 
central star like the Sun. Planets at a < 0.72 AU are Type II planets flHamano et al.ll2013l) . 
which are in steam planet state in the beginning and lose their water before cool down. Thus, 
they never have liquid water on their surface. On the other hand, planets at a > 0.72 AU 
and with Min, < 0.05 Mqc may be land planets from the beginning. 
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Fig. 7.— The evolution of planets with Mini = 1-0 Mqc for central stars like the Sun. None 
of these planets become land planets. They have too much water to start with. Note that 
some of the planets that orbit older stars spend considerable time as steam planets. 
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Fig. 8.— Like Figure [7] for initially drier planets with Mini = 0.1 M^c- Aqua planets between 
0.72 and 0.82 AU can transform into land planets, and as land planets they maintain liquid 
water on their surfaces for an additional 2 Gyrs. Thereafter, such planets receive the solar 
radiation of the threshold of the runaway greenhouse for a land planet and become steam 
planets. 
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Fig. 9.— The inner edge of the CHZ, dehned as the planet being in the HZ for at least 
4.6 Gyr, considering the evolution from an aqua planet to a land planet. The region of the 
type-I CHZ corresponds to the traditional CHZ for an aqua planet, while the region of the 
type-H CHZ is a region where initially aqua planets become land planets. 
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Fig. 10.— The effect of the critical vapor amount for the runaway greenhouse Mcv on the 
maximum initial amount of water, Mini to evolve from the aqua planet mode to the land 
planet mode. The maximum water mass for a land planet, which is the other transition 
condition, is hxed at Mmi = 0.05 Mqc. The result for (solid curve) the case where M^ is 
equivalent to a column of 3 m of precipitable water is the same as that shown in Figure |6l 
Larger critical amount of vapor facilitates the evolution to land planets. However, if is 
less than Im, it is very difficult for an aqua planet to evolve into a land planet. 
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Fig. 11.— The effect of the maximum water mass for a land planet, Mmi, on the maximum 
initial amount of water, Mini to evolve from the aqua planet mode to the land planet mode. 
Here Mcv, the other transition condition, is hxed to 3 m depth in global average. The result 
for the case of M^i = 0.05 Moc was shown in Figure [6l Changing M^i has a smaller effect 
on the upper limit of initial water mass to evolve from an aqua planet to a land planet than 
changing The amount of water for a land planet from the beginning is corresponded 

with the maximum water amount for a land planet. 





